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he violation of local uncertainty relations is a valuable tool for detecting entanglement, especially in 
I jmulti-dimensional systems. The orbital angular momentum of light provides such a multi-dimensional 
-—system. We study quantum correlations for the conjugate variables of orbital angular momentum 
^nd angular position. We determine an experimentally testable criterion for the demonstration of 
^ Tan angular version of the EPR paradox. For the interpretation of future experimental results from 
our proposed setup, we include a model for the indeterminacies inherent to the angular position 
' bieasurement. For this measurement angular apertures are used to determine the probability density 
^tjf the angle. We show that for a class of aperture functions a demonstration of an angular EPR 
C _ ^ paradox, according to our criterion, is to be expected. 

(y-jl, Introduction 

o , 

"^jExperiments on the orbital angular momentum [0AM] of light confirmed re- 
Qpently an uncertainty principle for angular position and angular momentum 
+IjiV/hereas for separable quantum states uncertainty principles limit the accu- 
racy for measurements of non-commuting observables, inseparable or entan- 
^gled states may apparently overcome these limits. This was first discussed by 
Qfiinstein, Podolsky and Rosen in their famous Gedankenexperiment 2 and 
^led to the formulation of the EPR paradox [2] . The implications of the EPR 
■^.^Jparadox have been tested mainly on optical systems, e.g. on the polarisa- 
''^ion of photons [UlEj, quadrature phase components jUEj or directly on the 
C^ptical version of EPR's original example, the linear momentum and linear 
position of photons |H]. The relation between 0AM and its conjugate variable, 
the angular position, is fundamentally different from these systems, because 
0AM is a discrete quantum observable of infinite dimension and the angular 
position is continuous and bounded. The entanglement for 0AM of photon 
pairs generated in parametric down conversion has been confirmed both ex- 
perimentally and theoretically jlUj . It is therefore interesting to examine 
the possibility to demonstrate an 'angular' EPR paradox for this pair of ob- 
servables, in particular as the necessary experimental techniques have already 
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been employed in recent work in this field [Tllllj. 

Apart from the fundamental interest in angular EPR correlations, the cri- 
terion for an EPR paradox for 0AM and angular position of light provides 
a tool to characterise entanglement for these observables. Using variances in 
violations of local uncertainty relations as entanglement criteria has received 
renewed interest. Formerly applied to continuous variables on a specific exam- 
ple [H], it has been shown that variances of special observables can be used 
to detect entanglement in finite-dimensional system This approach has 
been generalised to arbitrary observables and led to entanglement criteria for 
finite dimensional systems which are known from the continuous variable 
regime |14| . 

A practical motivation for studying entanglement of optical 0AM arises 
from the possible advantages of 0AM in quantum information processes (15] 
and quantum communication, as cryptographic schemes ^Hl could profit from 
an enlarged basis of states. 

1.1 Orbital angular momentum 

0AM of light is connected with the azimuthal phase structure of light beams: 
each photon in a beam with phase dependence exp{im(j)) carries an 0AM of 
mh \n. 18, 19i| . The 0AM number m G Z can take on any integer number 
which leads to a discrete, infinite dimensional quantum system. The conju- 
gate variable is the angular or azimuthal position (pQ E [9,2Tr + 6), which we 
choose to lie within the 27r radian interval starting at angle 9. The associated 
uncertainty relation is then PUIIT] 



where AL^ = hAm and P{9) is the angular probability density at the bound- 
ary of the chosen interval. The topology of the basis sets is reflected in their 
Fourier relation; in contrast to the linear case a discrete Fourier transform 
allows us to change from 0AM representation to angle representation: 
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Here, ip{(f)) is the wavefunction in the angle representation and Cm the 0AM 
probability amplitude. The 27r radian interval is commonly chosen to be 



Angular EPR paradox 



3 



[— tTjTt), and we set 9 = — vr from here on. 

It has been experimentally demonstrated that 0AM of light is conserved un- 
der parametric down conversion [5]. In theoretical studies, it has been pointed 
out that the conservation of 0AM is related to the phase matching condition 
for parametric down conversion processes ^Oj. Entanglement in 0AM and 
azimuthal position is a consequence of this phase matching The conser- 
vation of transverse momentum requires that the two-photon wavefunction for 
the signal (1) and idler (2) mode, for a plane wave pump, has to be of the 
form 5{ki^x + k2^)5{ki^y + k2,y)- Using a simplified approach one can argue that 
the transverse spatial correlations in the far field originate from the momen- 
tum conservation under parametric down conversion. Identifying transverse 
momentum components in the near field with spatial coordinates in the far 
field allows us to write the spatial dependence of the wavefunction in position 
representation as 

5{xi + X2)5{yi + y2) = —5{pi - P2)52t,{4>i - <;/'2 - vr), (4) 

Pi 

where p and (/> are the radial and azimuthal coordinates and 52-k is the 2tt 
periodic delta function. From this result we can expect the azimuthal angles 
of the photons in the far field to obey 0i = (/>2 + vr, so that the signal and idler 
photons appear on opposite sides of their respective cones. The correlation in 
angular momentum follows on writing Eq. Q in terms of its angular Fourier 
components: 

1 oo 

6{xi+X2)5{yi + y2) = —5{pi- P2)— (-1)™ exp(im0i) exp(-im02), 

' m=—oo 

(5) 

which is an entangled superposition of states with zero total 0AM. A more 
detailed analysis which considers a specific parametric down conversion pro- 
cess shows that more complicated dependences of the wavefunction on the 
azimuthal angles are also possible j21j . 

1.2 EPR paradoxes 

The EPR paradox describes the apparent violation of the uncertainty principle 
resulting from measurements on correlated, spatially separated systems. The 
original EPR argument considers correlations that are strong enough to predict 
or infer with certainty the value of the observables in one subsystem from 
measurements on the other, separated, subsystem without disturbing in any 
way the first subsystem. The ability to predict with certainty the value of an 
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observable defines, according to EPR, an element of reality. However, non- 
commuting observables, cannot have a simultaneous reality, an expression of 
which is the uncertainty principle [^. The tension between local elements of 
reality and quantum complementarity leads to the paradox. 

The original EPR Gedankenexperiment considers an idealised situation. The 
quantum state given by EPR on the example of the position and momentum is 
— in the modern language of entanglement - a maximally entangled state |23j . 
and the measurement of the observables is assumed to be infinitely precise. For 
0AM and angular position this idealised setting would require a parametric 
down conversion process which creates an entangled photon pair, perfectly 
correlated in 0AM and angular position. An errorless measurement of the 
0AM on the signal photon could then be used to infer the 0AM of the idler 
photon, and an errorless measurement of the azimuthal angle on the signal 
photon would allow us to predict precisely the angle of the idler photon. As 
these measurements on the signal photon 'do not disturb the idler photon in 
any way', the predictions would constitute simultaneous elements of reality 
for the 0AM and the azimuthal angle of the idler photon. We stress that the 
possibility to predict observables of the idler photon with certainty depends 
on the ability to measure the observables on the signal photon without error. 

In particular for a continuous observable a measurement with infinite pre- 
cision cannot be realised experimentally. A typical experimental setup would 
allow us to determine whether a continuous variable falls into a previously 
specified range. To analyse the possibility of demonstrating an EPR paradox 
experimentally a more realistic situation has to be studied. This requires the 
consideration of non-maximal correlations and of measurements with finite 
precision leading to an error in inferring one observable from a measurement 
on the other subsystem. The size of this error determines whether the EPR 
paradox can be demonstrated in the considered experimental setup |24j . 



2 Formulation of the paradox 

The inclusion of experimental indeterminacies requires a reformulation of the 
EPR paradox. In contrast to the original EPR setup 0, we consider condi- 
tional measurements on both subsystems. In reference to the idealised setup 
in section (|1.2|) . we look at the variance of the 0AM or the azimuthal position 
in the idler beam given a specifically set outcome of a measurement in the sig- 
nal beam. For the 0AM this condition will be the measurement of a single m 
value, whereas for the azimuthal position the measurement will be in a range 
of angles. For the 0AM we denote the conditional variance with var[m2|W'i], 
i.e. the variance of m2 in the idler under the condition that a measurement 
on the signal photon yields an 0AM of mi. The condition for the azimuthal 
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position to fall in a range of angles will be treated as a probability density 
Pi{4'1t'Ti)- The functional dependence of the probability density is given by 
Pi{(j)i) and the variable ri is used to indicate the orientation. For a symmetric 
Pi(0i), Tl would be the central angle. 

The conditional variance for the angle can thus be written as 
var[(^2|-Fi(</'i; "T"!)]- For an actual experiment the error caused by non-maximum 
correlations is thus included in measured quantities and the theoretical mod- 
elling, therefore, concentrates on the description of the error in the measure- 
ment. Schemes to measure 0AM have been theoretically studied and experi- 
mentally realised j25pilj . We are therefore mostly concerned with a description 
of the angle measurement. In the following we give the criterion for the exper- 
imental demonstration of an angular EPR paradox as one main result of this 
work. 

2.1 Criterion for an angular EPR paradox 

In the previous section (|1.2() we pointed out the importance of inference in EPR 
type arguments. To distinguish a measured quantity from an inferred one we 
label the first with an index m and the latter with an index i. The measured 
quantities are the conditional probability for the 0AM -P[m2|m'i]m and the 
conditional probability density for the azimuthal angle i'[</'2|-Pi(</'i; Ti)]m. The 
paradox becomes now apparent if one assumes local realism. From the mea- 
sured variance i'[</'2|-Pi('/'i; Ti)]m a minimum variance min var[m2|Pi((/>i; Ti)]i 
can be derived which is still in accordance with the uncertainty relation [cf. Eq. 
(^]. This quantity can be compared to the measured variance var[m2|?n.i]m 
as the conditioning measurement on the first subsystem 'does not have an 
instantaneous influence on the second subsystem'. An angular EPR paradox 
would then be demonstrated if 

(var[m2|mi]m)mi < (minvar[m2|Pi(0i; ri)]i)^^. (6) 

Within the simplified reasoning in section 11.11 the correlations in 0AM are 
uniform and the correlations in angle isotropic. This cannot be assumed a priori 
for an experimental test of the criterion Eq. © . We are therefore considering 
averaged conditional variances, which take into account that the correlations 
may vary for different values of mi and for different orientations of the angle ri. 
In the quantity (var[ 

77X2 |^i]m )mi the Conditional 0AM variance var[?7i2|^i]m is 
averaged over the condition mi. To find the average (minvar[m2|-Pi((/'i; ''"i)]i)ri 
the minimum conditional variance min var[m2|Pi((/)i; Ti)]i is integrated over ti 
on a 27r radian interval with the probability density P{(f)i]Ti). Formally the 
averaging procedure eliminates the dependence of the criterion on particular 
values of mi and ti, but not the dependence on the functional form of 
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The dependence on originates from the way in which the condition in 

Eq. @ is measured and is not directly connected to the correlation in the 
azimuthal angle. 

With the usual definition of the variance in terms of probabilities, we can 
reformulate the paradox statement: 



The inferred conditional probability P[m2\Pi{(pi] ti)]\ will be calculated via a 
Fourier transform from P[(j)2\Pi{(j)i; Ti)]m and is given by the modulus square 
of the conditional probability amplitudes c[m2\Pi{(j)i;Ti)]\. An experimental 
result which obeys the given inequality [Eq. would constitute a demon- 
stration of an angular EPR paradox. 

2.2 Angle measurement scheme 

The measurement scheme to determine the azimuthal position models exper- 
imental techniques employed in recent work J4. A photon is said to have a 
specific angular probability density if it is detected after passing an angular 
aperture corresponding to this probability density [see Fig. With the help 
of spatial light modulators smooth aperture functions can be realised. 

A narrow aperture can be used to measure the probability density of the 
azimuthal position. The aperture could then be rotated, i.e. the central angle 
could be varied over a 27r radian range. Detecting the number of photons 
passing the narrow aperture as a function of the central angle yields eventually 
a measure for the probability distribution. 

2.3 Conditional variances 

The angular apertures will also be used to set the condition in the signal beam. 
A fixed aperture can be inserted in the signal beam, and a rotatable aperture 
can be used to measure the conditional probability density in the idler beam 





(7) 




[see Fig.E]. 
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Figure 1. (a) Measurement scheme for the azimuthal position. A photon is said to have a 
particulaj: probability density for the aaimuthal angle if it is detected after passing an aperture 
corresponding to this probability density. Experimentally these apertures may be shaped using a 
spatial light modulator, (b) Aperture functions and their associated azimuthal probability densities 
for a transmitted photon, shown for a rectangular aperture and an aperture in form of a truncated 

Gaussian. 




-It It -Jt It 

♦ ♦ 



Figure 2. Schematic representation of the conditional measiurement. A fixed angular aperture in 
the signal beam sets the condition and a rotatable aperture in the idler beam can be used to 
measure a conditional probability density for the azimuthal position. The difference in the aperture 
functions is chosen here for the purposes of illustration only. 
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We will denote the outcome of such a measurement by P[(j)2\Mi{Ti)]m, 
i.e. the probability density to detect an idler photon with angle (p2 under the 
condition that the entangled photon in the signal beam passes the aperture 
Ml oriented at ri. As the aperture can be described by an aperture function, 
which we assume translates into a probability density Pi {(pi ; ri ) , we will write 
synonymously -P[02|Pi(0i; ''"i)]m- The importance of the orientation angle ti 
has been stressed in the formulation of the paradox [cf. section [2.1] . To sim- 
plify the notation we will not write the explicit dependence of Pi on ti from 
here on. 

From the conditional probability density a class of wavefunctions in the 
angle representation 'ip[(j)2\Pi{4'i)] is derived, 



i;[4>2\Pi{<Pi)] = VP[(t>2\PiiMUewiMh))- (8) 

Here, the phase ia((^2) is undetermined, as the measured probability densi- 
ties only give the modulus square of the wavefunction. The wavefunction is 
then transformed into a conditional 0AM probability amplitude via a Fourier 
transform [cf. Eq. 

1 



c[m2\Pi{(l)i)]\ = —7= d4>2 exp {i'm(p2 ) \/P [02 1 Pi (^i )] m exp(ia((?:)2 ) ) . (9) 
v27r J_7r 

From the conditional probability amplitudes we can calculate the conditional 
variance var[m2|P((?!)i)]i by taking the sum over all 7712 values 



min var[m2|P((/>i)]i = min 



E^,c[m2|Pi((/ >i)]P/^ 
c[m2|Pi( 



12 



Em. c[m2\Pl{<t>l)]fl 

V c[m2|Pi(<Ai)]P 



(10) 

This is the inferred minimum variance which can be compared to the measured 
quantity var[m2|7Tii]m- The phase ia{4>2) will be determined by the minimiza- 
tion of the conditional variance var[m2|Pi(<?5<i)]i as detailed in the following 
section. 



2.4 Minimization of the conditional variance 

Calculating the conditional wavefunction from the conditional probability 
leaves the phase ia(</)2) undetermined. We find that if the variance is calculated 
for a{4>2) = then the minimum variance minvar[m2|Pi((/'i)]i is obtained. To 
show this we assume the conditional wavefunction to be of the form 



ip[4>2\Pii<Pi)] = Aexp[ia(02)], 



(11) 
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where A = ^\ip[(j)2\Pi{4'i)]\'^ is a positive, real function, which is periodic 
in <j)2- Applying the orbital angular momentum operator to the wavefunction 
yields 

= -\hA'{(f)) exp[iQ(02)] + ha{(i)2) exp[ia((?!)2)], (12) 

where the primes denote derivatives with respect to (j)2- Using the periodicity 
of A the variance of L can be evaluated to 



varL, = h^ j #2^(02)^" (<^2) + j #2(a'(</'2))'P(02) - 

= #2a(</'2)P(02)) , (13) 

where we used the fact that A^{<f)2) = P[(t>2\Pi{4'i)\m- The first integral is 
the variance of L for a{(j)2) = 0, while the second and third integral are the 
variance of a': 

waiLz = [varLj;]a=o + ^^vara'. (14) 

Therefore we obtain the minimum variance if the wavefunction ip[(j)2\P {(j^i)] is 
real and positive, so that a = 0. In the following we will consider only this 
case. 



3 Proposed experimental scheme 

In order to measure var[m2|mi]m we have to examine the signal photon for 
the particular 0AM mi and determine the 0AM m2 of the idler photon. A 
schematic representation of an experimental setup to achieve this is shown in 
Fig. (121) . A spatial light modulator is used to produce a hologram which changes 
the 0AM in the signal photon by the chosen value —mi to zero |25| [T] . Only 
beams with mi = have on-axis intensity and can thus be detected behind 
a pinhole. The measurement for the idler photon has to distinguish between 
different values of m2. A sorting scheme which is able to determine the 0AM 
has been experimentally implemented TTl. A coincidence measurement would 
then yield var[m2|mi]m- 

Experimentally it will not be possible to measure the conditional proba- 
bility density for a single angle, instead a suitable aperture can be used to 
test for a range of angles, as described in section [T^ Analogously to the con- 
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pump 




Figure 3. Experimental scheme for measuring var[77i2|mi]m- In the signal a hologram is used to 
single out a particular value of mi as a condition. The OAM distribution of the idler is determined 
with help of an interferometric OAM sorter 1111 . Only the first stage of the sorter is shown here, 
additional stages are added where indicated by arrows. Each stage doubles the possible outcomes 
and therefore the number of detectors. Eventually the signals from all detectors are transmitted to 

a coincidence counter. 



ditioning aperture, the measured quantity is P[M2(T2)|Pi((/)i)]m) where the 
aperture M2 is centred at a particular angle (/)2 = T2. The aperture can be 
described by a probabiHty density P2i4'2',T2), where T2 indicates the central 
angle. The measured probabihty can thus be written as i'[^'2('/'2; ''"2)|-Pi(</'i)]m- 
For a very narrow aperture M2 this measurement will give a good estimate of 

P[r2|Pl(0l)]m 

P[r2|Pi(0i)]m ~ i^[i^2(</'2;r2)|Pi((/'i)]m for suitable P2(</'2). (15) 

The error made in this approximation is then given by the variance of 02 for 
the probability density P2. 

The scheme to measure P[P2{4'2',T2)\Pi{4'i)]m is basically shown in Fig. ((JJ. 
With the help of a spatial light modulator [SLM] any chosen Pi{(j)i) can be 
set as a condition for the signal. For the idler the aim is to determine (/>2 as 
exactly as possible. To achieve this an SLM could be programmed for a narrow 
angular aperture P2{4>2]T2) centred at (j)2 = T2, which would then be varied 
over the 27r radian interval. Eventually this would lead to P[P2{4'2)\Pi{(t^i)\m 
which does not only depend on the condition Pi{(f)i) but also on the chosen 
analysing aperture P2{(t>2)- This is the experimentally measurable estimate 
of the quantity i'[<?!'2|-Pi(</'i)]m used in the formulation of the angular EPR 
paradox [cf. section H2.1() ]. Obviously there are experimental limitations: a 
narrow aperture would transmit only little intensity, which would reduce the 
detection rate in the conditional measurement. Also SLMs have a finite size and 
resolution, which limits the ability to distinguish between similar apertures. 



Angular EPR paradox 



11 



The interpretation of experimental data would require the inclusion of im- 
perfect correlations originating from the parametric down conversion process 
and the influence of the optical elements, in particular the angular apertures. 
In the following we will discuss these aspects briefly. 

3.1 Parametric down conversion 

The possibility to create photon pairs entangled in 0AM or angular position 
is based on the conservation of 0AM under parametric down conversion |H1 
llUj . which holds for thin down conversion crystals and in the paraxial limit. 
The conservation leads to a perfect anti-correlation in the 0AM indices so 
that for a given 0AM index nip in the pump, signal and idler obey rup = 
mi + m2- In a recently reported down conversion experiment j^, the spatial 
correlations of photons entangled in orbital angular momentum have been 
studied. In this particular experiment signal and idler cone overlap completely 
and the spatial correlations are such that for a fixed detector position in the 
signal, the coincidence pattern in the idler shows two distinct spots equally 
separated in angle from the position exactly opposite the signal detector on the 
phase matching ring. The vertex of the separation angle is on the pump axis. 
In our work we are concerned with a different angle: the azimuthal position of 
a photon in a beam is measured from the beam axis, i.e. in a down conversion 
experiment from the signal and idler axis respectively. For non degenerate 
down conversion crystals these two angles are not identical. The question if 
they are compatible for the degenerate case is a very interesting one, but it 
will not be examined in the scope of this work. Therefore, for our theoretical 
modelling, we make use of the simplified reasoning given in section (|1.H) and 
therefore have an angle correlation — (/)2 = vr in the far field [cf. Eq. (jl}]. 

A recent study on the EPR paradox for linear optical momentum and po- 
sition [HI included a theoretical prediction of the conditional probabilities. 
Although this is certainly an interesting additional information for interpret- 
ing experimental data, the analysis here does not rely on it, as we would like 
to end with a criterion that can be applied directly to experimental data. The 
effect of imperfect correlations is therefore included in the measured quantities 
var[m2|mi]m and P[P2(02)| A(<Ai)]m- 

3.2 Angular apertures using spatial light modulators 

The advantage of spatial light modulators [SLMs] is that, within the spa- 
tial resolution of the SLMs, angular apertures may be smooth functions of 
the angle (j). Rectangular functions, which represent 'cake-slice' apertures, 
would result in singular derivatives and hence in an infinite inferred vari- 
ance var[m2|-Pi(0i)]i t27j. However, this analysis does not take into account 
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that optical diffraction will have a smoothing effect on the angle probabil- 
ity distribution. To study the influence of even very small smoothing effects 
we consider in the following theoretical modelling a class of continuously dif- 
ferentiable aperture functions which asymptotically approximate rectangular 
functions. On the other hand, aperture functions which differ only slightly may 
be mapped to the same aperture in the SLM, because of the limited resolution. 
These aspects have to be discussed more closely in conjunction with specific 
experimental implementations. 



4 Theoretical modelling 

To give a quantitative result we model the measurement of P[P2{4>2)\Pii<Pi)] 
under the two assumptions that the photon pair is perfectly correlated in angu- 
lar position and that the angle probability distribution behind the aperture is 
exactly given by the function describing the aperture. Under these assumptions 
the conditional probability P[P2(</'2)|-Pi(</'i)] is given by the overlap integral of 
the two probability densities, since the detection probability for an analyzing 
aperture centred at (j)2 = ^2 is given by 

PdU<p2=r2)= f I \ij{4>l,<t>2)?Pl{ct>l)P2{<p2:T2)d<Pid4>2. (16) 

J —TV J —TV 

Using the assumption of perfect correlation the probability density will be 
sharply peaked for (pi — (j)2 = and in that sense we may use a 27r-periodic 
(5-function to approximate |V'((/'i, 02)P 

|V^(0i,</.2)|2«52.(0i-</,2-7r). (17) 

The idler probability density P2{4'2 \ T2), may be written as a function centred 
at (/)2 = but shifted by T2, which yields ^2(02 — T2). Using these results the 
detection probability can be rewritten as 

/TT 
i^l(0l)i^2(01+Vr-T2)d,/.l. (18) 

For the calculation of this integral it is important to use the periodicity of the 
probability densities if the argument lies outside the 27r radian interval. The 
detection probability -Pdtc(</'2 = T2) is the measured conditional probability 
P[P2{4>2', T2)\Pi{(pi)]m- Varying T2 over the 2tt radian interval yields the condi- 
tional probability P[P2{<l)2)\Pi{4>i)]m- In an experimental test of our criterion 
[Eq. this quantity would be known from measurements. 



Angular EPR paradox 



13 



4.1 Aperture functions 

We model the conditional angle measurement for different aperture functions. 
As experiments to validate the uncertainty relation used apertures which 
can be described by truncated Gaussians, we calculate the effects of these 
truncated Gaussian apertures and also of a set of truncated super Gaus- 
sians. The latter allows us to interpolate between Gaussian and rectangular 
apertures. As mentioned in section H3.2() rectangular apertures, which could 
be made from a solid, absorbent material, will lead to an infinite inferred 
variance aperture functions discussed in this section are symmetric, 

i.e. Pj{(pj) = Pj{—(j)j) for j = 1,2. The overlap integral in Eq. (|18l )can there- 
fore be turned into the convolution of the two probability densities 

Pi{-ct>i)P2{(l)i + vr - r2)#i = [Pi * P2](vr - T2). (19) 

-TT 

For simplicity we are modelling the quantity P[(t)2\Pi{4'i)\m with the convo- 
lution [Pi * P'2\{(t)2) as for the symmetric and periodic aperture functions the 
shift by TT and the sign of T2 is not relevant. 



4.1.1 Rectangular functions. The probability functions describing rectan- 
gular apertures can be given in terms of Heaviside step functions H[(j)): 

P^{(t>j) = —H{4>j + ^)Hi-(^j + ^) for j = 1,2. (20) 

According to Eq. (|19j) the conditional probability is given by the convolution 
of the two probability densities 

Ai - A2 Ihl < A2, 
Ai -(^.2 A2 < < Ai, (21) 
Ai < \(t>2\ 

where we introduced the notation Ai = {wi + W2)/2 and A2 = {wi — W2)/2 
for wi > W2. As the convolution is commutative [Pi * P2 = P2 * Pi] we can 
interchange wi and W2 for the case W2 < wi. The conditional wavefunction is 
then given by the positive square root. Graphs of the rectangular apertures 
and the resulting conditional probability density and wavefunction are shown 
in Fig (01). 

The Fourier integral [Eq. ©] can be calculated analytically by using the 



[Pl*P2](02) = 
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Figure 4. In a) the probability densities representing the apertures in the signal and idler beam 
are plotted, for a width in the idler and signal of wi = 1/Att and W2 = l/647r respectively. The 
resulting conditional probability density and wavefunction are shown in b). The wavefunction is 

taken as the modulus of the probability density as a purely real wavefunction minimizes the OAM 



variance. 



Fresnel sine and cosine integrals ^2 and C2 j28j : 



c[m2 1 Pi (</>!)] 



1 



VA1-A2 



'MlI^C2i\m2\{A,-A2)) 



cosdmal Al) 



52(|m2|(Ai-A2)) 



(22) 



The distribution of conditional probability amplitudes is given in Fig. (jSJa) for 
the analytical calculation and a numerical integration. The numerical results 
are shown to give an estimate of the accuracy of our integration. This serves 
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as a reference for aperture functions, where we have not found an analytical 
solution of the Fourier integral. The conditional variance is shown in Fig. (Ob) 
over a maximum 0AM index m2 at which the sum in Eq. 1)10(1 is truncated. 
The numerical results differ from the analytical points at higher truncation 
indices because of numerical effects in sampling the 27r radian interval. Since 
the Fresnel integrals tend to 1/2 for large arguments, the conditional ampli- 
tudes c[m2|Pi(</>i)] vary with |m2|~^ for large m2, which leads to a logarithmic 
increase with higher truncation indices and thus to a divergent conditional 
variance. The reason for this behaviour is founded in the singular derivative 
of the conditional wavefunction at (j)2 = Ai. 



4.1.2 Truncated Gaussians. Aperture functions in form of truncated Gaus- 
sians have been used in the experiment studying the angular uncertainty prin- 
ciple ill. In this case the probability densities Pi and P2 are 

PM) = -r / ^ ^^P ( - f— ) ) ^ t-^'^)'-?' = (23) 

V^'u;jerf(7r/-u;j) \ ywj J J 

The graph for Gaussian apertures is shown in Fig. ^a) for the same widths 
wi and W2 as in the rectangular case. The resulting conditional probability 
density and wavefunction are plotted in Fig. ^h). 

The conditional probability density P[(j)2\P{4>i)] can be analytically calcu- 
lated using the convolution theorem for the Fourier transform 



1 



-\/27r J-TT 



1)2 exp(im202)[Pi * P2](02) 



1 



1 



27r erf(7r/?i;i)erf(7r/w2) 



X exp 
xRe 



'm2{wi + W2 



vr .m2Wi 
err 1 

wi 2 



(24) 



Re 



erf ( ^ - i^] 

\w2 2 ; 



An approximate analytical expression for the probability amplitudes can be 
obtained by treating the Gaussians as extended which results in setting the 
error functions to unity. This allows us to calculate the c[m2|-Pi((/'i)] from Eq. 



\ vr / \erf(7r/t(;i)erf(7r/tt>2) /V ^ J 

(25) 
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b) 


numerical 256 + ^ 




numerical 512 x 




analytical □ 




logarithmic fit 

+ 




X 

+ X □ 




fi ^ 








20 40 60 80 100 120 140 160 
truncation 



Figure 5. a) Numerically and analytically calculated OAM distribution for rectangular apertures. 
The analytical solution is exact, but to give an estimate of the accuracy of our integration method, 
the numerical results are also shown. In b) the conditional variance is plotted. The numerical 
results are given for the two different sampling sizes 256 and 512. The deviation of the numerical 
256 results for higher truncation indices is caused by numerical effects. For the analytical points the 
logarithmic increase in the variance can be seen. 



For the chosen widths the agreement with the numerical solution is excellent 
as can be seen in Fig. ((7|a). This is, however, a singular case and for different 
parameter settings the approximation differs more from the numerical solution. 

The resulting conditional variance is shown for numerical and analytical cal- 
culations. For the truncated Gaussians the variance converges. The difference 
between the numerical results and the analytical solution is relatively small, 
because the contributions from larger values of m2 are very small. 
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Figure 6. In a) the probability densities representing the apertures in the signal and idler beam 
are plotted, for a width in the idler and signal of wi = l/4n and W2 = l/647r respectively. The 
resulting conditional probability density and wavefunction are shown in 6). The wavefunction is 
taken as the modulus of the probability density as a purely real wavefunction minimizes the OAM 

variance. 

4.1.3 Truncated super Gaussians. With truncated super Gaussian [TSG] 

apertures we can gradually go from the rectangular case to the Gaussian. This 
is achieved by an additional parameter 7 in our definition of the TSG: 



= X — Jti — TT^iTT exp I 1 for G [-tt, tt), j = 1, 2, 




(26) 

where r(-) is the complete Gamma function and r(-; •) the incomplete Gamma 
function. For values of 7 > 1 we have a kurtosis smaller than 3, i.e. probability 
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Figure 7. a) Numerically and analytically calculated OAM distribution for truncated Gaussians 
apertures. The analytical solution is approximated by neglecting the finite interval and treating the 

truncated Gaussians as infinite. This allows us to calculate the probability amplitudes using the 
convolution theorem and the fact that square roots of Gaussians are again Gaussians. The resulting 
conditional variance can be seen in b). For the Gaussian apertures we have a constant conditional 

variance. 



densities, which are less peaked than Gaussians. The effect of this parameter 
on the probabihty density can be seen in Figs. (|Hla) and ((Hlb). 

For the TSG all calculations have been done numerically and one can see 
m Fig (Pa) that, even for high values of 7, the probability distribution, 
|c[m2|Pi((/>i)]p, differs substantially from the rectangular case. For the vari- 
ance, the dependency on the truncation index m2 and 7 is shown in Fig. (jUJb). 
For small values of 7 the variance does not change over the range of trunca- 
tion indices, for higher values the variance converges more slowly or not at 
all within the given range. For the highest value of 7 = 80 the effect is more 
pronounced, but far from the logarithmic increase in the rectangular case. 
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Figure 9. a) Numerically calculated 0AM distribution for TSG apertures. The distribution is 
plotted for two different values of 7. In b) the dependency of the conditional variance on the 
truncation index m2 and 7 is shown. One can see that for small 7, i.e. 7 < 5 the variance converges 
and stays independent of the truncation index. For higher 7 the variance converges more slowly, 
but there is no logarithmic divergence. 

4.2 Discussion 

The theoretical modelhng of a conditional angle measurement leads to a con- 
ditional variance significantly different from zero. However, the final result 
depends not only on the aperture in the signal, which sets the condition, but 
also on the aperture used to determine the angle. The angular EPR criterion 
Eq. ® takes the orientation of Pi {(pi ; ti ) into account by averaging over the 
the orientation angle ri. From the simplified approach in section [TTT] we could 
assume isotropic correlation for the azimuthal angle, which is expressed in the 
52tt function in Eq. Q. Within this model the correlation are therefore the 
same for every orientation ti and an averaging is not necessary. The depen- 
dence of P[P2{(j)2)\Pii4'i)] on the aperture functions P2 and Pi reflects the 
way in which the azimuthal angle is measured. 
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From the presented results, it is clear that the broadness of the conditional 
angle probability density and the 0AM probability distribution is mostly de- 
termined by Pi{4>i)- The influence of the analyzing aperture -P2(02) is certainly 
most controllable for the rectangular case. But the analysis of TSGs showed 
that there are values of 7 for which the influence of P2{4>2) is comparable 
to the rectangular case without showing the divergent variance. Using these 
apertures to study the influence of the analyzing aperture on the conditional 
probability distribution could lead to a more detailed model for the measure- 
ment process and thus to a more complete theoretical prediction for the final 
conditional 0AM variance. On the other hand, for TSGs the question remains 
how well the aperture can be programmed in an SLM, in particular as for high 
values of 7 the features which distinguish it from a rectangular aperture are 
on a very small scale. This however might be an interesting aspect in a more 
detailed analysis of rectangular apertures, which includes optical diffraction, 
as even small deviations from a perfect rectangular form would give a finite 
conditional variance. 



5 Conclusion 

In this work we have discussed the possibility to demonstrate an angular EPR 
paradox for the conjugate variables of orbital angular momentum [0AM] and 
angle. The paradox is about the apparent violation of an uncertainty rela- 
tion for incompatible observables measured on correlated, spatially separated 
subsystems. By using photon pairs entangled in 0AM and in angle these sub- 
systems can be realised in an optical experiment. We have found a testable 
criterion for an angular EPR paradox, which takes experimental indetermi- 
nacies into account. For that we have reformulated the EPR paradox using 
conditional variances, i.e. variances of observables from one subsystem given 
a preset outcome on the other subsystem. 

To investigate the feasibility of an experimental demonstration, we have 
modeled the measurement process under the assumption of perfect angle cor- 
relation. Also, angular apertures to set the condition or to determine the angle 
were assumed to impose their probability characteristics exactly on transmit- 
ted photons. Under these assumptions we have studied different classes of 
aperture functions for the final conditional 0AM variance. Rectangular func- 
tions lead to a divergent conditional 0AM variance, which does not set any 
lower bound for the correlations in 0AM for our criterion. Truncated Gaus- 
sians result in a quickly converging variance, but the implementation of the 
apertures with a smooth grayscale transition on an SLM will be less accu- 
rate than for rectangular apertures with their sharp edge. Truncated super 
Gaussians, which can be varied from the rectangular case to the truncated 
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Gaussians, provide an aperture which leads to a convergent variance and still 
has a controllable influence of the analysing apertures in the idler. The condi- 
tional variances obtained from the theoretical modelling show that an angular 
EPR paradox can be demonstrated. Given the current state of experiments 
we expect an implementation of our criterion to be able to demonstrate the 
EPR paradox for 0AM and azimuthal position. 
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